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PACS 42 . 50 . Ex - Optical implementations of quantum information processing and transfer 
PACS 42.50.Ct - Quantum description of interaction of light and matter; related experiments 
PACS 42 . 50 . -p - Quantum optics 

Abstract. - We propose a method that enables efficient storage and retrieval of a photonic excita- 
tion stored in an ensemble quantum memory consisting of Lambda-type absorbers with non-zero 
Stokes shift. We show that this can be used to implement a multimode quantum memory stor- 
ing multiple frequency-encoded qubits in a single ensemble, and allowing their selective retrieval. 
The read-out scheme applies to memory setups based on both electromagnetically-induced trans- 
parency and stimulated Raman scattering, and spatially separates the output signal field from the 
control fields. 



The distribution of entangled states over large distances 
is a basic requirement for realizing quantum networks and 
numerous quantum communication protocols [1,2]. It is 
also of importance for fundamental experimental tests of 
quantum theory [3] . As recently pointed out [4] current en- 
tanglement distribution protocols using quantum repeater 
stations between two parties [5,6] will only work over large 
distances if significantly improved quantum memories be- 
come available. Ideally these memories should also be 
able to store multiple qubits and allow their selective re- 
trieval [7]. Here we propose a scheme to achieve these 
goals for photonic qubits stored in a quantum memory 
(QM) consisting of an ensemble of three level A-type ab- 
sorbers as shown in Fig. 1(a). 

Our scheme is based on the approach [8-14] where the 
qubit is encoded in a signal light pulse, and stored in the 
ensemble via a classical control field [see Fig. 1(a)] propa- 
gating colinearly with the signal in the positive z-direction 
[as defined in Fig. 1(b)]. Storage of the signal produces a 
highly-asymmetric spin wave in the medium [15, 16]. Af- 
ter some time the pulse can be retrieved by applying an- 
other control field. If this field also propagates in the pos- 
itive z-direction, the spin wave of the medium has a low 
overlap with the read-out mode, resulting in a low mem- 
ory efficiency unless one resorts to a much larger coupling 
strength for read-out. Instead one can reverse the read- 
out control field, which maximizes this overlap. However 
for non-zero Stokes shift W13 [see Fig. 1(a)] the memory 
efficiency suffers due to momentum conservation. It is de- 
sirable to have LO13 ^ for independent addressability of 
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Fig. 1: Scheme for efficient storage and phasematched 
spatially-resolved read-out in an ensemble QM. (a) Absorber 
level configuration, and (b) schematic experimental setup. A 
signal pulse (frequency lo 3 , wavevector k s ) is stored by a con- 
trol field (frequency uj c , wavevector k c ). The configuration is 
such that initial state | 1) is energetically higher than storage 
state I 3) . The control field is angled so that the resulting spin 
wave has zero z-momentum, allowing efficient retrieval of the 
signal in a different direction to the read-out control field. 



the transitions | 1) <-» | m) and | 3) <-> |m), as selection 
rules are rarely stringent enough to ensure this. Further- 
more, in colinear quantum memories one also has to spec- 
trally separate the retrieved signal field from the stronger 
control field. For a weak signal this filtering is difficult, 
and another method for resolving the fields, such as ori- 
enting the control field slightly off-axis [17, 18], would be 
advantageous. 

The aim of this paper is to overcome these limita- 
tions in the read-out. Our proposal applies to QMs 
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based on both stimulated Raman scattering [15,19,20] and 
electromagnetically-induced transparency (EIT) [8,11,13, 
14] . We show that using our method, a multimode QM can 
be realized whereby frequency-encoded qubits are stored 
in and retrieved selectively from a single ensemble. This 
would allow manipulation of stored qubits, for example in 
an optical lattice QM [21], without needing to use mul- 
tiple ensembles for the different logical states. The num- 
ber of modes that can be stored in this way is limited 
experimentally by angular resolution and the achievable 
frequency of signal and control fields. In a conventional 
single-mode memory the signal and control pulses are col- 
inear and the levels are chosen so that state 1 1) is energet- 
ically lower than | 3). In this case W13 7^ 0, and the stored 
spin wave has a momentum Ak = k s — k c (we call this 
a phase mismatch) , rendering backwards read-out of the 
signal inefficient. To overcome this (i.e. to phasematch 
the system) we choose the level configuration shown in 
Fig. 1(a), and orient the read-in control field at an an- 
gle of 9 C = cos _1 (|k s |/|k c |) to the signal propagation di- 
rection [see Fig. 1(b)]. This eliminates the longitudinal 
phase mismatch and allows efficient retrieval of the signal 
field. Furthermore the output signal field direction, de- 
termined by momentum conservation, is spatially distinct 
from the read-out control field. The scheme allows high- 
fidelity retrieval of signal pulses with duration T provided 
that L <C Tc, and Tuiis ^> 1, with L the length of the en- 
semble and c the speed of light. In the following we detail 
the application of this method to a general QM, including 
a discussion of the constraints that arise and the storage 
times that can be achieved. We then show how our scheme 
could be used to implement a multimode Raman QM con- 
sisting of a single ensemble, and derive an expression for 
the maximum number of modes it is possible to store. 

We consider the setup shown in Fig. 1. The classical 
control field for read-in is oriented in the x-z plane at 
an angle of 9 C to the signal. If 6 C is small, the walk- 
off between the two fields can be neglected. The con- 
trol field envelope can then be represented at time t and 
longitudinal position z by the slowly varying Rabi fre- 
quency fi(r), where r = t — z/c. The frequency of the 
transition | 1) <-> | m) is Lj\ m . We introduce dimension- 
less coordinates £(z) = z/L for the longitudinal position; 
e(r) = u(t)/w(T), with uj(t) = /J" |Q(r')| 2 di"' the in- 
tegrated Rabi frequency, for the time; and p = (X, Y), 
where X = x^J uj\ m /2Lc (analogously for Y) for the 
transverse position. We define the slowly-varying oper- 
ators a(e,C,p) = y/uj(T)/CA{r, r)e-^^ z ) /fi( r ) for the 
signal field and /3(e,(,p) = ^/tjCB{v, r) e - ix(r '^ for 
the spin wave, where A(r, r) and B(r,r) are the slowly- 
varying signal and spin-wave amplitudes respectively, and 
x(t, z) = [w(r) + |K| 2 ,z]/r, which describes a Stark shift 
due to the control field and a modification of the signal 
group velocity. We decompose the phase mismatch as 
Ak = Ak^ + Afc 2 z, with z the unit vector oriented along 
the z-axis. The coupling of the signal field to the ensemble 
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Fig. 2: Demonstration of the increase in read-out efficiency for 
an ensemble QM with reverse read-out. Retrieval probability 
v is plotted against the read-in and read-out couplings Ci n 
and Cout (defined in text). In (a) we use colinear signal and 
control fields. The phase mismatch introduced by non-zero 
LJ13 degrades the efficiency, with v — 0.23 for Cm = C ou t = 2, 
which is sufficient for storage with near-unit efficiency [15]. In 
(b) the control field is angled appropriately, and d n = C ou t = 2 
gives v = 0.95, which is also much higher than the efficiency 
achieved using colinear read-in and colinear forward readout 
{v = 0.52) [15]. 

is |k| 2 = dy/L, with d the resonant optical depth [16] and 
r = A — i7, where 7 accounts for dephasing and loss of 
the optical polarization. 

We assume that either or A are much larger than the 
signal bandwidth, the maximum control Rabi frequency, 
and 7, and adiabatically eliminate the excited state | m). 
Within this limit our theory describes the interaction for 
arbitrary signal and control pulse shapes in both the EIT 
and Raman regimes. The absorbers are all initially opti- 
cally pumped into state 1 1), and the population of state 
I 3) is assumed to remain negligible for each absorber. 
In the slowly-varying and paraxial approximations the 
Maxwell-Bloch equations are 

d e (3 = -Ce- ipX a, (1) 
(5+i9 C )« = Ce^A (2) 

where q = |k s |c/cJi TO , p = |k e | ain(6 c )(uJi m /2Lc)~ 1 / 2 , and 
the coupling C = \n\^J Luj(T) /\T\. Note that these equa- 
tions are similar to those in [15], except we also include 
the transverse Laplacian V 2 = d\ + d\ [22]. 

The dynamics of Eqs. (1) and (2) decomposes into a 
linear mapping between signal and spin wave modes. The 
optimally-efficient mapping excites an asymmetric spin 
wave Pi approximated by exponential decay in £. The 
symmetry of the interaction [15, 16] dictates that the opti- 
mal spin wave for forwards retrieval, (3q, is the mirror im- 
age in £ of this, i.e. /3q(C) = /Ml — C)i where a superscript 
r identifies a quantity associated with the read-out. The 
efficiency of forward readout then depends on the overlap 
of (3i with /3q, which is generally low. Switching the prop- 
agation direction of the control field, to readout in the 
backward direction, flips f3 v around, so that its overlap 
with /?i is perfect. However, due to the residual momen- 
tum Ak of the spin- wave coherence after storage, we then 
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have ffi (C, p) = Pi ((, p) exp [up (£, p)], where 
¥>(C,p) = {^^(Ak x + Aki).p+(Afc, + A^) 

xiC + x[r,iC]-x r [o,i(i-C)]}- 

The linear ^-dependence of the real terms in ip reduces 
the overlap of the stored spin wave and the output mode, 
leading to a reduction in the memory efficiency. This is 
shown in Fig. 2(a) for a cesium QM, where levels 1 1) and 
| 3} are the F = 3 and F = 4 hyperfine levels of the 6S , 1 / 2 
state respectively (i.e. | 1) is energetically lower than | 3)). 
Our strategy is to hnd a configuration for which 3?{<p} is 
independent of £. To do this we choose 1 1) to be ener- 
getically higher than | 3) and set 9 C = 6% = cos _1 (k s /k c ), 
where the angle of the read-out control field 9 T C is as de- 
fined in Fig. 1(b). This phascmatches the memory in the 
Raman and EIT configurations - far from resonance the 
X functions in (p are real but negligibly small, and on res- 
onance they are imaginary, so do not contribute to the 
mismatch. Figure 2(b) shows that this results in a high- 
fidelity read-out in 133 Cs, where now | 1) and | 3) are the 
F = 4 and F = 3 hyperfine levels respectively. Note that 
a transverse phase mismatch does not necessarily degrade 
the readout efficiency; any residual Ak^ simply modifies 
the angle at which the retrieved signal field is emitted. 
However, having 9 r c ^ 9 C requires a modification of the 
frequency of the control field to ensure that momentum 
and energy are conserved. 

We have assumed that the walk-off along the z-axis be- 
tween the signal and control fields is negligible, i.e. small 
compared to the pulse durations, so that 9 C <C \jTcjL. 
A more stringent condition is that the transverse walk- 
off of the beams is small compared to their beam waists, 
which requires 9 C <C VA/L, with A the cross-sectional 
area of the control pulse. As well as setting the angle 
9 C , the Stokes shift imposes an upper limit to the control 
field bandwidth, since there should be no control photons 
at the signal frequency. Diffraction limits the length of 
the ensemble according to the condition A ~ AL, where A 
is the control wavelength. These conditions can be com- 
bined into Tc > L > VA, with the maximum storage 
bandwidth given by TCJ13 ^> 1. The upper bound on L 
is relaxed if the control field focussing is loosened, but the 
control pulse energy should then be increased. However, in 
general, the above considerations imply that the number 
density required for high-efficiency (>90%) phasematched 
operation of a A-cnscmble QM is given approximately by 
n = 6T~ 2 , with 6 ~ lnr 3 s 2 , in both EIT and Raman 
configurations. 

In a 133 Cs QM the 65*1/2 state is split by the hyper- 
fine interaction with the nuclear spin (/ = 7/2), to give 
the F = 3 (our |3)), and F = 4 states (our 1 1)), with 
Wi3 = 9.2 GHz . The 6P 3/ / 2 state (D? line), which lies 
351.7 THz (852 nm) above state 1 1), would be used for the 
excited state | m). If A = 10 GHz, then the phasematch- 
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Fig. 3: Storage and retrieval of two signal modes with different 
frequencies: (a) the stored spin wave intensity and (b) the 
retrieved signal field intensity I a . Both are Fourier-transformed 
in the ^-coordinate to illustrate the different spectral compo- 
nents. The two signal modes require control fields angled at 
9 C — ±3°, thus creating two components to the spin wave at 
either side of k x = 0. Read-out of the two modes is performed 
simultaneously with a control field, consisting of two frequency 
components, angled at 9 T C — 3°. One output photon propagates 
along the z-axis, and the other propagates at an angle of 6° 
to the z-axis. Alternatively two separate control fields, one af- 
ter the other, could be used to retrieve the modes at different 
times. 

ing angle is 9 C ~ 0.5°. This would enable high-efficiency 
storage and retrieval of a single photon with T = 250 ps 
in an ensemble with L ~ 2 cm. This bandwidth requires a 
number density of around 10 19 m~ 3 , which can be achieved 
in a cesium vapor at a temperature of T e ~ 360 K. For 
a thermal vapor we must consider motional dephasing of 
the memory, as the phase introduced by the read-in con- 
trol field will average out [23,24]. The stored spin wave 
oscillates at the frequency clk-lsin^ ~ c|k c |# c , and so 
the entanglement fidelity of the memory [25] is given by 
T ~ exp[— \k c \ 2 9 2 t 2 s kBT e /M], with t s the storage time of 
the memory, ks Boltzmann's constant, and M the atomic 
mass. For the above parameters the maximum storage 
time with T > 0.9 is t s ~ 200 ns, i.e. three orders of 
magnitude larger than the signal pulse duration. Longer 
storage times could be achieved by using a storage unit 
consisting of an ultra-cold gas, but this requires a larger 
sample than has currently been achieved experimentally 
Using the off-axis QM scheme described above, it is pos- 
sible to store multiple signal fields of different frequen- 
cies in a single ensemble, and retrieve them selectively, 
as demonstrated in Fig. 3. Hence one can implement a 
multimodc QM [7] in which frequency-encoded qubits are 
stored, and retrieved either one after the other, or at the 
same time but spatially separated. To sec this we con- 
sider the same setup as in Fig. 1, but now the signal field 
consists of two components, labeled 1 and 2, of frequency 
lu Si i and lo s _2 respectively. Signal component j (j = 1, 2) is 
Raman resonant with control field j , which has frequency 
w C j, detuning Aj, and is angled at 9 c j to the signal for 
phasematching. Both signal components can be stored in 
the ensemble, and the resulting spin wave has two compo- 
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ncnts separated in k x (where k x is the momentum com- 
ponent in the x-direction) , due to the different transverse 
momenta of the control fields [see Fig. 3(a)]. Backwards 
read-out of each mode can be performed using a control 
pulse with the appropriate orientation and frequency, in 
analogy with the single mode case. These control fields 
could be applied at different times to read out the modes 
one after the other, or the angles can be changed so that 
the modes are emitted in different directions [as shown in 
Fig. 3(b)]. We ensure that |Ai - A 2 | > (5i + 62), with 
Sj the bandwidth of signal component j, so that signal 
field j cannot initially interact with control field i (i 7^ j). 
This condition requires that max[Ai,A2] is large, i.e. in 
the Raman regime. However, once a material excitation is 
created by the absorption of signal j, this excitation may 
interact with control field i. 

To calculate the form of the stored spin wave and the 
retrieved signal field, we again assume that the angles rel- 
ative to the signal field of both control fields arc small. 
Control field j has wavevector in the x-dircction with 
rescaled magnitude pj. Wc assume that control field j 
has Rabi frequency £lj(r) = f2j/(r), with f2j the maxi- 
mum Rabi frequency, /(r) an envelope function normal- 
ized so that Uj{T) = \ flj\ 2 T, and the integrated Rabi fre- 
quency uij(r) = /J" |ri J (r')| 2 dr'. Also due to the large 
number of absorbers in our ensemble, the probability of a 
single absorber being excited by both signal components 
is assumed to be negligible. The coupling constant is now 
C m = kVLW, where W = £^1 ^( r )/l r jf- The di- 
mensionlcss operator for signal component j is written as 
0j (e, C,p) = y/T/C^e-^^Ajjr r)//( r), and the spin 
wave operator is now b(e, (, p) = \J LjC m e^ lXb{ - T ^B(v, r), 
where Aj(r,t) is the slowly- varying amplitude of signal 
field j, and the spin- wave momentum mismatch is now 
Ak = ZW13/C. The complex exponents are Xj( T > z ) = 
YLiUi(j)/T t + MVr,, and X b(t,z) = EL<*W/r? 
respectively. The Maxwell-Bloch equations are thus 

ii^l + i9 cK = C mCj e-^ x+ ^b, (3) 
d e b = -C* m [e^ x+ Wc iai + e^ x+ V R ^c 2 a 2 ] , (4) 

with Cj = y/T/WiSlj/Tj), Cj = ^T/W{n*/Tj) and 
Rj = Tj^/W /\k\VL; subscripts indicate the signal com- 
ponents to which the quantities refer. The solution for the 
two-component signal field, which can be achieved in the 
same way as for the single-field case, is the convolution 
of the integral kernels for the single-field interaction [15] 
with terms that arise due to the mixing of the two com- 
ponents. This mixing occurs, for example, when during 
storage signal component j interacts with both control 
fields by a second-order process, and could potentially de- 
grade the read-in. However, our calculations have shown 
that the mixing only becomes significant when the read-in 
coupling is C m > 15 - much larger than the coupling suf- 
ficient for signal storage with near- unit efficiency (C = 2). 
For C m = 2 we see that the storage process for the two- 
component signal field is almost identical to the results 



obtained when the two signal components are assumed to 
interact independently of each other. The read-out pro- 
cess also obeys Eqs. (3) and (4), and a similar lack of 
mixing ensures that the modes can be addressed indepen- 
dently Note that if the cu's have a relative phase of n 
between them, then the two absorption processes interfere 
destructively, and the storage efficiency is low. For effec- 
tive storage of the signal components this phase must be 
zero, which can be ensured either by changing the sign of 
the Aj's, or by phase-shifting the Slj's. 

Since one can in principle have arbitrarily-far blue- 
detuned fields (avoiding any higher levels), the main lim- 
itations on the number of modes that can be stored are 
experimental - the achievable optical frequencies and the 
possible angular resolution AO. The maximum control 
field angle allowed is given by 6* max = s/A/L and the num- 
ber of different control fields allowed is [0 max /A0\ . The 
smaller angles considered here correspond to the fields be- 
ing far blue-detuned, which may be outside the optical 
regime. Taking this into account, the number of modes 
that can be stored is given by n m = [{(9 max — [2^13/ (ct>i TO + 
A max )] 1/ ' 2 }/A6'J , with A max the largest allowed blue de- 
tuning from I m). For the 133 Cs QM discussed previously, 
taking the angular resolution to be AO = \/\/~A, setting 
A = 10~ 7 m 2 , and restricting A max so that no other tran- 
sitions are excited, this corresponds to n m ~ 100 modes. 
As with the single-mode case the phase introduced by con- 
trol field j is LD c ,jO c .j / c = uJxs/ (c6 C) j). Hence the limita- 
tion on storage time of the memory due to dephasing, as 
calculated earlier for a single mode, is given by the en- 
tanglement fidelity of the mode that requires the smallest 
control field angle for phascmatching. 

In summary, we have demonstrated a QM scheme using 
an ensemble of three-level A-type absorbers that allows ef- 
ficient storage and phase-matched spatially-resolved read- 
out of the signal field. We use this method to demon- 
strate how a multimode QM could be implemented using 
a single atomic ensemble. Multimode memories will en- 
able the storage of multiple qubits in one memory setup, 
significantly decreasing the resources required for entan- 
glement distribution protocols. Furthermore, storing mul- 
tiple qubits in a single ensemble opens up new possibilities 
for scalable entanglement manipulation. 
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